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Abstract. Motivated by recent experiments on multi-component membranes, the 
growth kinetics of domains on vesicles is theoretically studied. It is known that the 
steady-state rate of coalescence cannot be obtained by taking the long-time limit of 
the coalescence rate when the membrane is regarded as an infinite two-dimensional 
(2D) system. The steady-state rate of coalescence is obtained by explicitly taking 
into account the spherical vesicle shape. Using the expression of the 2D diffusion 
coefficient obtained in the limit of small domain size, an analytical expression for the 
domain growth kinetics is obtained when the circular shape is always maintained. For 
large domains, the growth kinetics is discussed by investigating the size dependence 
of the coalescence rate using the expression for the diffusion coefficient of arbitrary 
domain size. 
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1. Introduction 

Lipid bilayer membranes can be regarded as two-dimensional (2D) systems embedded in 
three-dimensional (3D) solvent. The membranes are coupled to solvent since the lipids 
composing the membrane interact with solvent surrounding it and the momentum can 
be exchanged between the membrane and the solvent [Il[2l[3j|Il[5l|6j[7J[8j|9j ITU] . 
In this sense, membranes can be regarded as quasi-2D systems. By the recent 
advances in experiments, domains formed by phase separation in multicomponent 
membranes are visualized, and the lipid spatial organization and its dynamics have been 
studied [HJ H2J EEl HU [15] . The phase separation kinetics and coarsening are influenced 
by the dimensionality, domain shapes and the hydrodynamics in the systems [16]. The 
phase separation in multicomponent lipid bilayer exhibits rich dynamics partly due 
to the momentum dissipation to the third- dimensions while the motion is confined to 
2D [7J El E]. 

In recent experiments, the growth kinetics of circular domains on ternary fluid 
vesicles has been observed by fluorescence microscopy [TH [T3J, [TJ] . In these experiments, 
liquid domains are formed in giant vesicles by phase separation into the liquid-ordered 
phase and the liquid-disordered phase on lowering the temperature from the one- 
phase region. Yanagisawa et al. found two different types of domain coalescence 
kinetics [H]. In one of the coalescence processes, the domains grew by collision and 
coalescing while keeping the circular shape until the large domains of the vesicle size 
appeared. This growth kinetics due to the diffusion-controlled coalescence (DCC) was 
described by a power-law. In the other coalescence process, the domain growth was 
suppressed by membrane-mediated repulsive inter-domain interactions. Recently, it 
was pointed out that the domain coalescence could be prevented by the membrane- 
mediated interactions between liquid domains associated with the deformations of the 
membrane [TTJ, [TEJ ElHl EH EE]. When the liquid domain size exceeds a critical value, 
the boundary line energy is reduced by budding at the expense of some bending energy. 
The significant slowing down of the domain growth can be observed for large budded 
domains [HI [21]. In the steady state, the domain patterns and the membrane shapes 
can be stabilized by the coupling between the local membrane curvature and the line 
tension pi [151 E2 E31 EI] . 

As briefly summarized above, the liquid domains coarsen under the influence of 
several competing processes. Even without budding, the observed growth kinetics was 
different from that obtained from the scaling hypothesis [131 E3]- Motivated by the 
experiments, we study theoretically the growth kinetics of domains on vesicles by DCC. 

The study of domain growth on vesicle surfaces is still limited compared to that 
in 3D |16j . According to the scaling hypothesis, the domain growth exponent due to 
DCC in 2D is 1/2 in contrast to that in 3D given by 1/3 [251 EH]. However, it should 
be noted that, in the scaling hypothesis, the coalescence process is not explicitly taken 
into account. Moreover, it is known that the steady-state rate of coalescence cannot be 
obtained for an infinitely large 2D system [271 EH] . In contrast to diffusion in 3D space, 



Growth kinetics of circular liquid domains 



3 



the boundary conditions are crucial to obtain the coalescence rate when the diffusion is 
restricted in 2D. This point is missing in the scaling hypothesis. 

In addition to the above argument on pure 2D systems, the coalescence of liquid 
domains can be influenced by the coupling between the membrane and the solvent. The 
diffusion of large domains is more influenced by the coupling compared to that of small 
domains. As the domain grows, the influence of the coupling between the membrane 
and the solvent increases. The domain growth kinetics has been studied by dissipative 
particle dynamics simulations and continuum simulations [TJ El El EH]. The simulation 
results suggest the slowing down of the domain growth by DCC due to the coupling [7]. 

In this paper, we study the growth of liquid domains immersed in a 2D membrane 
by using an analytical theory which goes beyond the scaling hypothesis. We note that 
the steady-state rate of coalescence can be obtained by taking into account explicitly the 
vesicle shape [30J, [31] . By using the diffusion coefficient of domains obtained by taking 
into account the coupling between the membrane and the solvent [H [2j [3j HJ El [10] , 
we show that it is independent of the domain size in the limit of small domain size. 
In such a case, it is known that the size distribution is described by the Smoluchowski 
theory of coalescence processes [32]. By further assuming that the circular shape of 
the liquid domains is always maintained, the time evolution equation of the mean 
domain size is obtained from the size distribution using the conservation of domain 
area upon coalescence. The mean domain growth is expressed by a single function for 
the whole time regime starting from the initial induction period of coalescence to the 
final asymptotic regime given by the power-law. When the domain size is large, we 
discuss the influence of the coupling between the membrane and the solvent on the 
domain growth by analyzing the size dependence of the coalescence rate. 

In Sec. [21 we present the known results obtained from the scaling hypothesis. In 
Sec. [31 the results of Smoluchowski theory in pure 3D and 2D infinite systems are 
reviewed. In Sec. HI the coalescence rate is obtained by taking into account the vesicle 
shape. In Sec. El the analytical expression representing the growth of mean domain 
size is obtained when the circular shape and the area of the liquid domains are kept 
before and after the coalescence. The size dependence of the the coalescence rate is 
investigated by using the analytical expression of the diffusion coefficient for the liquid 
domain of the arbitrary size in Sec. El In Sec. the theoretical results are discussed in 
relation to those obtained by the recent experiments [14j . 

2. Scaling theory 

In this section, we summarize the results obtained from the scaling theory. The scaling 
theory is the simplest way to derive the power-law growth of domain size. Obviously, one 
cannot obtain both the transient growth leading to the asymptotic power-law kinetics 
and the magnitude of the power-law growth. 

The scaling theory is based on the hypothesis that the mean domain radius (a(t)) 
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Figure 1. Schematic picture liquid domains embedded in a fluid vesicle (left). The 
radius of liquid domains is a while the vesicle radius is R. Liquid domains undergo 
Brownian motion on the vesicle surface. Shown on the right is the geometrical 
parametrization. 9 denotes the azimuthal angle and a is the angle between domain 
centers at collision, a represents the encounter distance. 

in d- dimension is related to time t by 

(a(t)) 2 ~ D d t, (1) 

where Da is the diffusion coefficient and subscript d denotes the dimensionality. In 
3D, the diffusion coefficient of the domain is inversely proportional to (a(t)) by the 
Stokes-Einstein relation, D% ~ l/(a(t)). By substituting this relation in the scaling 
relation, we obtain (a(t)) ~ t 1 ^ 3 [16| |25| 126] . For 2D, we have D 2 ~ k^T '/r/ by the 
dimensional argument, where fee is the Boltzmann constant, T the temperature, r] the 
2D membrane viscosity surrounding the domain. By using the above expression of D2, 
the scaling hypothesis leads to (a(t)) ~ t 1 / 2 [T6J |25l [26]. It should be pointed out 
that a logarithmically time-dependent diffusion coefficient was derived in pure 2D by 
the hydrodynamic theory [2J. When the membrane couples to the solvent, a constant 
diffusion coefficient was derived, and D2 depends logarithmically on the domain size 
in the weak coupling limit [U |2]. The logarithmic size dependence and the numerical 
factors are ignored in the dimensional argument. 

We also note that the time given by {a(t)) 2 /Dd is not equal to the coalescence time 
of the domain with size (a(t)). The coalescence time should be given by the mean first 
time that domains coalesce by diffusion from random initial distribution. Below, we 
show that the latter time is different from the one given by (a{t)) 2 / D a- 

The scaling hypothesis merely relates a single length scale given by the mean domain 
size at time t with the diffusion coefficient as shown in Eq. ([1]) and it should be justified. 
In the subsequent sections, we point out that the results of the scaling hypothesis cannot 
be obtained for the domain growth by the diffusion-coalescence in 2D infinite systems. 
Then, we show that the scaling hypothesis is consistent with the domain growth kinetics 
by diffusion-coalescence when the available diffusion area is finite in 2D. The whole 
kinetics including transient growth and the final power-law growth will be obtained by 
taking into account the diffusion-coalescent process and the vesicle shape explicitly. 
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3. Smoluchowski theory in infinite systems 

In the Smoluchowski theory, domain motion is assumed to be independent of the others 
and coalescence between a pair of domains is considered. When both domains can move, 
it is difficult to solve the problem analytically. We assume that one of the domains is 
fixed and the other diffuses with the mutual diffusion coefficient D d which is the sum of 
the diffusion coefficients of two spherical domains of equal size [33J. Coalescence takes 
place immediately when two domains come in contact at the encounter distance which 
is the sum of the radii of the two domains. In the theoretical formulation, the spatial 
distribution of domains satisfies the diffusion equation and the spatial domain density 
should vanish at the encounter distance. The domain size increases immediately after 
the coalescence and the spatial distribution of the new domain size should be zero at 
the encounter distance for the increased domain size. 

The density of mobile domain around the immobile domain satisfies the diffusion 
equation 

^ t p d (r } t) = D d W 2 Pd (r } t) } (2) 

where r is the distance to the center of the immobile domain. We assume random initial 
condition given by, 

p d (r,t = 0) = 1. (3) 

The boundary condition applied at the encounter distance a is given by 

p d ( r = a,t) = 0. (4) 

We should set another boundary condition such that the density at infinite separation 
is unity, i.e., 

lim p d (r,t) = I. (5) 

r— >oo 

3.1. Three- dimensions (3D) 

Before investigating the domain growth on the 2D spherical surface, we shall briefly 
present the known results for 3D infinite systems and show that the corresponding 
results do not hold for 2D cases. 

For 3D infinite systems, the density profile which satisfies both of the boundary 
conditions is p%{r) — 1 — (a/r) [33]. The mean coalescence rate is given by the inward 
flow of domains across the surface at a 

k 3 = Ana'D, = 4*<tD s . (6) 

In a 3D fluid of viscosity r] s , the diffusion coefficient is inversely proportional to the size 
of the diffusing object, D 3 = k-^T /37rrj s a. By using this Stokes-Einstein relation, the 
mean coalescence rate can be expressed by k 3 = 4k B T/3r] s . Notice that a/2 is the radius 
of the spherical domain [32] . 
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3.2. Two- dimensions (2D) 

By applying the boundary condition at a, we obtain 

p 2 (r)=C\n(a/r), (7) 

where C is a constant determined from the other boundary condition at r — > oo. 
However, it is impossible to fix C because p 2 (r) diverges for r — > oo as long as C is 
finite. Unlike 3D infinite systems, the density cannot satisfy both boundary conditions, 
Eqs. @ and (jSJ), simultaneously. 



4. Smoluchowski theory in spherical surface 



The difficulty mentioned in the previous section for 2D can be overcome if the available 
diffusion area is finite. The density profile depends crucially on the shape of the 2D 
region. In the experiments [H], DCC of circular domains was observed on the vesicle 
surfaces. In principle, the coalescence rate can be obtained by modifying the method 
shown in the previous section applicable to the spherical region but the calculation is 
rather complicated. In this paper, we employ an alternative method. 

In the method, we investigate the life time of the density of mobile domains survived 
from collision to the immobile domain. In a confined region, the decay of the density 
can be well approximated by a single exponential. The time constant of the exponential 
decay can be reasonably obtained from the mean first-passage time of a mobile domain to 
the periphery of the immobile domain by assuming uniform distribution for the starting 
point. For the coalescence, the mean first-passage time is the mean coalescence time 
corresponding to the encounter time between two domains. When the initial position 
of the mobile reactant is z = cos 6 (see Fig. [T]for the geometry), the mean coalescence 
time t(z) satisfies the following equation (see the Appendix for the derivation) [30] 
D2_d_ 
R 2 dz 

The boundary conditions are 

/ dr(z 



:i-z 



d~z T{z) 



-1. 



(8) 



T(Z 



— cos a 



0. 



dz 



0. 



(9) 



z=l 



where a is the angle between domain centers at collision (see Fig. [T]). Notice that a 
simple geometric argument gives sin(a/2) = a/{2R) and cos a = 1 — a 2 /(2R 2 ) [30J. 
Equations (JS} and (jUJ) can be easily solved, and we obtain 

'2R 2 



t(z) 



D 2 



1 + z) 



a' 



By averaging over the random initial distribution, we obtain [30 



dz t(z) 



2 - a 2 /(2R 2 



D 2 



In 



2R 



a 



(10) 



1- (ct/2R) 2 

The probability that the mobile domain has not reached the immobile domain up 
to time t is given by exp (— t/r av ) when domains are initially distributed uniformly 
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in the spherical surface. The same probability can be expressed by using the bulk 
bimolecular rate as exp (— k^t/A), where the surface area available for diffusion is 
A = 2nR 2 [2— a 2 / (2R 2 )}. Comparing these expressions, we find that the bulk bimolecular 
rate can be calculated from the mean coalescence time by k 2 = A/t^ [301 EU EI]- 

With the use of the bimolecular rate, the mean coalescence rate is obtained 
as [301 EH El] 

4vr[l - {a/R) 2 ]D 2 {a) 



2\a{R/a) - 1 + {a/R) 2 
4nD 2 (a) 



(12) 
(13) 



2\n(R/a) - 1' 

where a <C R is used to obtain the second equality. The mutual diffusion coefficient is 
now expressed by D 2 (a) since it depends on the domain radius a as shown below. 

Unfortunately, the full size-dependence of the diffusion coefficient is not known 
for spherical vesicles. However, the analytical expression is known for a circular liquid 
domain which has the same viscosity as the outside of the domain for 2D flat membranes. 
In this case, the mutual diffusion coefficient was obtained by De Koker as [3] 

2k B T r Jf( z ) 

D 2 (a) = / dz^— r, 14 

nr] Jo z 2 (z + ua) 

when the two circular domains have the same radius a. In the above, J\{z) is the 
Bessel function of the first kind, rj is the 2D membrane viscosity, v = 2r] s /r] with rj s 
being the viscosity of the outer fluids. Here we have assumed that the viscosities of the 
liquid inside and outside the vesicle are the same. The analytical expression after the 
integration can be expressed using Meijer G-functions [HIE]. In the case of va <C 1, the 
above expression reduces to [6j [5] 

— Jn(— )-7+T 

where 7 = 0.5772 • • • is Euler's constant. Equation ( TT5|) is slightly larger than the mutual 
diffusion coefficient of the Saffman-Delbriick (SD) theory derived for solid domains under 
the condition of va <C 1. Strictly speaking, the SD result was also obtained for 2D flat 
membranes. Recently, it was shown that the SD result was applicable for spherical 
vesicles when a < R < \ jv [36J . 

By substituting Eq. (jT5l) into Eq. (|T3|) . we obtain 

kWT , . 

k 2 w — , 16 
V 

when a/R 1. The coalescence rate is independent of a in this limit. Then k 2 t 
approximately represents the area explored by a diffusive object of radius a during time 
t [3U ES]. Hence k 2 t is given by D 2 (a)t times the effective collision cross-section. The 
size independence is the result of the two opposing effects; with increasing the domain 
size, the diffusion coefficient decreases while the effective collision cross-section given by 
47r/[21n( J R/a) - 1] in Eq. (CGI increases. 



D 2 (a) 



(15) 
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Strictly speaking, the coalescence rate depends on the domain size if the coalescence 
occurs between domains having different sizes. However, such a size dependence is small 
for 2D because of the weak logarithmic dependence in Eqs. (Tl3|) and ( fl5|) . and will be 
ignored hereafter. In fact, the size dependence was not taken into account even to study 
the coalescence processes in 3D [32] . 

Equation ffl2|) is valid even if the condition a/R <C 1 is not satisfied. The full 
size-dependence of the coalescence rate will be studied by substituting Eq. ffT4|) into 
Eq. ( fl2l as we shall discuss in Sec. El 



5. Growth kinetics of liquid domains 

Now we consider the formation of m-fold domains from the initial domains with the 
same size. When the coalescence rate is independent of the domain size as discussed 
above, the number density of m-fold domain n m {t) at time t is given by [32] 

n m {t) (t/r) m_1 



m (l + t/r)' 



where 



k d n x 

and n\ is the initial number density of the primary domains, and kd the coalescence rate 
in d- dimensional space. 

In general, if the aggregates of large sizes are in the solid state, there are many 
possible shapes of aggregates. The shape characterization and the size distribution 
have been studied in the conventional theory of Brownian coagulation. Here, we study 
Brownian coalescence of liquid domains under the condition that the circular shape is 
always maintained. For comparison, we also consider Brownian coalescence of liquid 
domains in 3D infinite systems, where the spherical shape is always maintained. Under 
this assumption, the time evolution of the domain size can be obtained analytically. 

By assuming that the domain shape is immediately restored after the coalescence, 
the radius of the m-fold domain a m is determined by the area (d = 2) or volume (d = 3) 
conservation relation, 

a d m = maf. (19) 

By averaging a m over the distribution of m-fold aggregate [see Eq. ifPTj) ]. the mean 
domain radius (a(t)) at time t is obtained as 

m. = £ m m*M = L Li . 1/d (J/L.) . (20) 

In the above, we have used the polylogarithm function defined by [37] 

oo k 

Li s (*) = £--- (21) 
k=i K 
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Figure 2. The average domain size {a(t))/a\ as a function of dimensionless time t/r 
when domains are confined in the spherical surface. The solid line is obtained from 
Eq. (|20|) and the dashed line represents the asymptotic behavior given by Eq. (|22|) . 



In the asymptotic limit of t/r S> 1, Eq. (120]) reduces to either 
(a(t)) ^ (t\ 1/2 



CLi 2 \T, 

for the coalescence in 2D spherical surfaces, or 



(22) 



W«»~rYiV±r (23) 



for that in 3D infinite systems. Here T(x) is the gamma function. The time evolution of 
(a(t))/ai is shown in Fig. [2] for the domain coalescence on the vesicle. The asymptotic 
time dependence is well approximated by Eq. (1221) when t/r 3> 1. We can see the 
induction period of coalescence when t/r < 1. The induction period is characterized by 
the inverse of the apparent coalescence rate given by the coalescence rate times the initial 
number density of the primary domains, Eq. ( Fl8|) . By taking into account explicitly the 
vesicle shape, the finite coalescence rate is obtained from the Smoluchowski theory and 
the results show the induction period before the asymptotic growth. 



6. Size dependence of the coalescence rate 

In the previous section, the growth kinetics of the circular liquid domain by diffusion 
coalescence is obtained when the coalescence rate is independent of the size of the liquid 
domain. The result for the domain growth on the vesicle is obtained in the limit of 
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Figure 3. k2r)/ (k-gT) as a function of the dimensionless domain size va. The solid 
line is obtained from Eq. (|T^1) using the expression of the diffusion coefficient for a 
circular liquid domain which has the same viscosity as the outside of the domain. 
Rv = 5, 10 and 50 from top to bottom. The dashed lines are obtained by substituting 
Eq. ([15]) into Eq. The dots are obtained by substituting Eq. into Eq. (TT2)) . 




« 1. Below, we study the size dependence of the coalescence rate using Eqs. (|T2|) 
and (TH1) . The diffusion coefficient Eq. ( IT^l) is obtained for a circular liquid domain of 
arbitrary size. In the limit of va <C 1, the domain size dependence of the diffusion 
coefficient is logarithmic as shown by Eq. ( fl5l) . The coalescence rate is given by k^T '/rj 
independent of the domain size because of the cancellation of the logarithmic domain 
size dependence both in Eqs. (|T2l) and ffT5]) . In the opposite limit of pa 3> 1, the diffusion 
coefficient is obtained as [51 E] 

* < 24 > 

which is inversely proportional to the domain radius a. 

The coalescence rate obtained by substituting Eq. (|T3|) into Eq. ( JT2|) is shown by 
the solid lines in Fig. [3j For relatively small vesicle, say vR = 5, the coalescence rate 
is almost independent of the domain size and is approximately given by k^T/i]. The 
situation is very different for large vesicles. When vR > 10, the size dependence of 
the coalescence rate can be ignored only when va < 0.1. In the same figure, we show 
the results obtained by substituting Eq. ( fl5l) into Eq. ( fl2|) using dashed lines. As long 
as the the domain size dependence of the diffusion coefficient can be approximated by 
Eq. (TI5T) . the size dependence of the coalescence rate is weak. 
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When vR = 5, the size dependence is weak even beyond the regime given by ua <C 1. 
However, when uR > 10 and ua > 0.1 the coalescence rate decreases by increasing the 
domain size. In the figure, we show the results obtained by substituting Eq. f l24l) into 
Eq. (fT2l) using dots. The coalescence rates with the full size-dependence of the diffusion 
coefficients [Eq. (fl4l) ] approach the results shown by dots as the size increases. The dots 
rapidly decrease by increasing the domain size because of the power-law dependence 
of the domain size in Eq. ff24|) . When vR > 10 and va > 0.1, the coalescence rate 
decreases by increasing the domain size because of the strong size dependence of the 
diffusion coefficient. If the coalescence rate decreases by increasing the size, the growth 
kinetics of the domain size should be suppressed. The growth law given by Eq. ( 1201) and 
the result shown in Fig. [2] can be applicable when vR = 5, but should be applicable 
only when ua <C 1 if the vesicle radius is large such as vR > 10. 

We also note that the coalescence rate increases by increasing the domain size if the 
domain size becomes close to the vesicle radius. This dependence originates from the 
denominator of Eq. ( 1121) and can be interpreted as the finite size effect of the domains 
confined in the spherical surface. 

The diffusion coefficient Eq. ( Tl4|) shows the logarithmic size dependence [Eq. ( fl5l) ] 
and the power-law size dependence [Eq. $M§] as ua is varied. The dimensionless quantity 
va characterizes the coupling between the embedding bulk fluid and the membrane for 
given size. The coupling is small when raC 1. The domain growth and its asymptotic 
limit given by Eqs. (120]) and ff22l . respectively, are obtained in the weak coupling limit. 
In the strong coupling limit ua ^> 1, the domain growth is suppressed as a result of the 
small coalescence rate compared to that in the weak coupling limit as shown in Fig. |2j 
The results are consistent with the simulation results that the growth law is suppressed 
by increasing the hydrodynamic coupling when the governing mechanism is DCC in 
2D [HE]. 

The above conclusion is not altered if we use the diffusion coefficient for the solid 
domains instead of Eq. (THI) . Recently, the diffusion coefficient of the solid domain 
in the 2D flat membrane is approximated by a closed-form empirical expression [10J. 
The interpolation formula of the diffusion coefficient reproduces the logarithmic size 
dependence obtained from SD theory when w « 1. When the interpolation formula 
is introduced, the coalescence rate Eq. ( f!2|) is slightly smaller but the overall size 
dependence is similar to that obtained by using Eq. ( fH|) . The coalescence rate decreases 
by the decrease of the diffusion coefficient of the solid domain compared to that of the 
liquid domain of the same size. In general, the diffusion coefficient of the solid domain 
is smaller than that of the liquid domain of the same size since the friction between the 
membrane and the solid edge is larger than that between the membrane and the liquid 
domain. 

As shown above, the coalescence rate is constant in time when domains are confined 
on a spherical vesicle. The finite coalescence rate constant can also be calculated by 
considering a circular flat sheet of radius L by setting the absorbing boundary condition 
at r = a and the reflecting boundary condition at r = L (a < L). When a circular 
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domain of radius a is placed at the center of a circular flat sheet, the coalescence rate 
is given by Eq. (1131) by replacing R with L [38]. Within this substitution, the domain 
size dependence of the coalescence rate is the same as that in a spherical vesicle. The 
diffusion coefficient of a domain placed at the center of a circular sheet of radius L was 
obtained by taking into account the viscosity of a circular flat membrane pp. However, 
the result was limited for a <^ L and obtained by ignoring the hydrodynamic coupling 
between the membrane and the solvent. 

7. Discussion and conclusion 

We discuss the relevant length scales in the recent experiment on domain growth kinetics 
by diffusion coalescence [H] . In the experiments, giant vesicles with a diameter of about 
20 jttm undergo phase separation at 30 °C after the temperature drop from the one-phase 
region (42-43 °C). The circular domains were observed when the size exceeded an optical 
resolution of the microscope (roughly 0.8 /mi). In one of the coalescence processes, the 
large domains of the vesicle size appeared within several minutes. The domain growth 
by collision and coalescence was observed. In the other coalescence process, the domain 
growth was suppressed for several 10 minutes. The former process can be theoretically 
studied by assuming DCC. 

The value of v — 2r] s /r) can be estimated by using the typical values of r] s = 10~ 3 
Pa-s, and r\ given by 0.1 Pa-s times the membrane thickness 5 nm as v = 4.0 x 10 6 
m -1 . We estimate vR ~ 40 by introducing the typical radius of the vesicles 10 /mi. 
When vR = 50, the coalescence rate is almost independent of the domain size as long 
as a < 25 nm obtained from the condition of ua < 0.1. The domain size is much smaller 
than the optical resolution of the microscope such as 0.8 /mi. Therefore, when the 
domain size grows and reaches to the optical resolution, the coalescence rate decreases 
with increasing the domain size and the domain growth can be suppressed compared 
to that given by Eq. ( 1221 . According to Fig. [3l the coalescence rate is almost constant 
over the wide range of the domain size when vR = 5. This value corresponds to the 
vesicle radius close to 1 /mi which may be the maximum vesicle radius to observe the 
power-law growth given by Eq. ff22l . 

It should be remembered that the full size- dependence of the diffusion coefficient for 
spherical vesicles is not known. In this paper, the size dependence of the coalescence rate 
has been discussed by substituting the known expression of the diffusion coefficient for 
2D flat membranes. The full size-dependence of the diffusion coefficient for the spherical 
vesicle is needed to further develop the coagulation theory for the large domains within 
the current limit of optical resolutions. 

In the original work by Yanagisawa et al, the best fitted exponent 2/3 was obtained 
for the power-law domain growth [TJ] . However, the mechanism which leads to this large 
exponent is not well-understood. An attractive interaction between domains seems to 
be present due to the hydrodynamic flow around domains, which would accelerate the 
domain growth [39J. 
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Power law growth of circular domains can be induced by transport of molecules 
from one domain to another through the medium [161 HOI SD H2] . The growth of large 
domains is associated with evaporation of small domains, which is known as Ostwald 
ripening. According to the Lifshitz-Slyozov- Wagner theory of Ostwald ripening, the 
power law exponent is 1/3 [10]. In this paper, we studied the domain growth by the 
DCC mechanism and did not consider the evaporation and condensation mechanism in 
2D. We just remark that the exponent 1/3 is independent of the dimensionality and 
holds also in 2D [HI EJ H2] . 

In conclusion, we have investigated DCC mechanism for growth kinetics of the 
liquid domains on the fluid vesicles. By applying the bimolecular reaction theory in 
the spherical surface and using the the 2D diffusion coefficient, the 2D coalescence 
rate is found to be independent of the liquid domain size if it is small enough. As 
a result, the domain size distribution is given by the classical Smoluchowski theory. 
When the circular shape is always maintained, we have obtained the mean domain 
size for the whole time range [Eq. (1201) ] by using the domain size distribution and the 
area conservation relation. In the asymptotic long time limit, we expect the power-law 
behavior with the exponent 1/2 [Eq. (1221)]. 

The domain growth kinetics has been derived under the condition that the 2D 
coalescence rate is independent of the domain size. The condition is investigated by using 
recently obtained analytical expression for the diffusion coefficient of arbitrary domain 
size. When the vesicle radius is small, the coalescence rate can be well-approximated as 
a constant over the wide range of the domain size. When the vesicle radius is large, the 
coalescence rate becomes independent of the domain size only in the limit of the small 
domain size. In general, the coalescence rate decreases by increasing the domain size 
up to a certain size where the finite size effect dominates. The results are discussed in 
relation to the recent experimental observations of DCC in vesicles. 
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Appendix A. The mean coalescence time 

In this Appendix, we briefly present the derivation of Eq. (jSJ) when the density satisfies 
Eqs. (J2I) (EJ) . The density can be expressed using the probability of finding a pair of 
domains at the relative position r at time t if their initial relative position was 7*i and 
was uniformly distributed 




(A.l) 
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We introduce the survival probability that the pair has not coalesced up to time t if 
their initial relative position was r 

w d (r,t) = J dr f p d (r f ,t\r,0). (A.2) 

Because p d (rf, 0|r, — t) satisfies the backward Kolmogorov equation and hence 
p d (r f , t\r, 0) = p d (r { , 0\r, -t), we obtain jl3l IM 135] 



^w d (r,t) = D d V 2 w d (r,t), (A.3) 
with the initial condition 

w d (r,t = 0) = 1, (A.4) 
and the boundary conditions 

w d (r = a,t) = 0, lim w d (r,t) = 1. (A. 5) 

r— >oo 

Since 1 — w d (r,t) is the probability that the pair coalesce at time t, the mean 
coalescence time r(r) is given by 

roo roo 

r(r)= dtt—[l-w d (r,t)] = dtw d (r,t). (A.6) 
Jo at Jo 

In general, the mean coalescence time is called the mean first-passage time. Integrating 

Eq. (1A.3j) over time and using Eq. (1A.4|) . we find that the mean coalescence time satisfies 

Eq. ([8]) when domains are confined in the spherical surface. Note that in Eq. (IA.3I) . 

the diffusion equation in the spherical surface can be written by assuming azimuthal 

symmetry such that 

where z = cos 6. 
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